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Abstract 

We calculate the nucleon-nucleon scattering amplitudes in the ^Sq and ^S^i-^I^i 
channels at next-to-next to leading order starting from a recently proposed non-relativistic 
chiral effective theory, which includes dibaryon fields as fundamental degrees of free- 
dom. We restrict ourselves to center of mass energies {E) smaller than the pion mass 
{tUt^)^ and further divide the calculation into two relative momentum {p) regions, a 
high energy one p ^ m^^ » 5m^, Snii being the dibaryon residual masses, and low 
energy one p < Srui. We first match to a lower energy effective theory in which we 
calculate the amplitudes in the high energy region. We further match this effective 
theory to the so called pionless effective theory in the low energy region, and carry 
out the calculations in the latter. Dimensional regularization and minimal subtraction 
scheme are used throughout. For ^Sq channel a good description of the phase shift 
data is obtained for E < 50MeV. For the ^Si Di channel, the ^Si phase shift data 
is only well described up to £" < 20MeV. 

PACS: 14.20.Pt, 13.75. Cs, 21.30.Fe, 21.45.Bc, 03.65.Nk . 
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1. Introduction 



Since the original suggestion by Weinberg [1] that the nuclear forces could be under- 
stood within the framework of effective field theories (EFT) there has been an enormous 
development of the subject (see [2, 3, 4, 5, 6, 7] for reviews). A key ingredient of the EFT 
formalism is that the cut-off dependence which is introduced in order to smooth out ultravi- 
olet (UV) singularities can be absorbed by suitable counterterms, and hence any dependence 
on physical scales much higher than the ones of the problem at hand can be encoded in a 
few (unknown) constants. In order to achieve this in a systematic manner counting rules are 
also necessary. 

In a recent paper [8] we proposed a chiral non-relativistic EFT which included two 
dibaryon fields as fundamental degrees of freedom. This EFT, which will be simply referred 
as NNEFT in this paper, is renormalizable and has simple counting rules when dimensional 
regularization (DR) and minimal subtraction (MS) scheme are used. The nucleon-nucleon 
scattering amplitudes in the ^Sq and ^Si channels were calculated at NLO and a good de- 
scription of data achieved in the — 50MeV energy range. We carry out here the calculation 
at N^LO in order to see if the good description of data persists and check the convergence 
of the EFT. This is mandatory in view of the fact that the so called KSW approach [9] also 
produced a good description of data at NLO, but turned out to have a bad convergence in 
the ^Si channel at N^LO [10, 11]. We will restrict ourselves to an energy range E such that 
E <C m^^ the pion mass, and p = ^/mNE ^ m^^ ttin being the nucleon mass. Pion fields 
can then be integrated out leading to an EFT, which was already described in [8], which 
we will call potential NNEFT (pNNEFT). For p - this EFT will already be suitable to 
carry out the calculations of the amplitudes. For p <^ ttIt^ however it will be convenient to 
integrate out nucleon fields with p ^ and use the so called pionless NNEFT (/NNEFT) 
[12, 13, 14]. All matching calculations will be done expanding the low energy or momentum 
scales in the integrals and using DR to regulate any possible IR divergence. Local field re- 
definitions which respect the counting will be used to get ride of redundant operators, rather 
than using the on-shell condition. 

We will organize the paper as follows. In sections 2 and 3 we introduce the NNEFT, 
and the pNNEFT Lagrangians respectively. In section 4 we match NNEFT to pNNEFT. In 
section 5 we calculate the nucleon-nucleon scattering amplitudes in pNNEFT. In section 6 
we match pNNEFT to /NNEFT, and calculate the nucleon-nucleon scattering amplitudes 
in the latter. Sections 7 and 8 are devoted to the comparison of our results with data in the 
^Sq and ^Si-^Di channels respectively. We close with a discussion and conclusions in section 
9. 
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2. The nucleon-nucleon chiral effective theory with 
dibaryon fields 



Our starting point is the effective theory for the Nb=2 sector of QCD for non-relativistic 
energies much smaUer than recently proposed in [8] . The distinct feature of this EFT is 
that in addition to the usual degrees of freedom for a NNEFT theory, namely nucleons and 
pions, two dibaryon fields, an isovector (D^) with quantum numbers ^Sq and an isoscalar (Dy) 
with quantum numbers ^Si are also included. Since ttiat ^ A^, a non-relativistic formulation 
of the nucleon fields is convenient [15]. Chiral symmetry, and its breaking due to the quark 
masses in QCD, constrain the possible interactions of the nucleons and dibaryon fields with 
the pions. The A^^ — sector is given by the chiral Lagrangian, which will only be needed 
at LO, 

-f 2 a a 

C^ = i^{Tr{d^U^d''U) + 2BoTr{MU^ + UM^)}, U = e'"^ (2.1) 
8 

Ad is the quark mass matrix, which we will take in the isospin limit, namely the average of 
the up and down quark masses rUq times the identity matrix. is defined by = 2Bomq. 
The Nb = 1 sector contains the pion- nucleon interactions , and will be needed at NLO [1], 

^ttN = + ^nN • (2-2) 

At LO we have, 

Ci'>=N^{.D,-,An-l)^£^)N (2.3) 

where ^ U, Uj^ = i [u'^.d^u^, = (9^ + ^[u^d/^u]), ir^ is the pion field, the isospin 
Pauli matrices, gA ^ 1.25 is the axial vector coupling constant of the nucleon, and /^^ ^ 
132MeV is the pion decay constant. This is the leading order Lagrangian for the pion- 
nucleon interactions {0{p)) for E ^ /2mN ^ m^- The NLO Lagrangian in this sector 
reads [2], 



A^ 

with x+ = 2Bq{u'^ Aiu^ + uA4^ u) . All parameters, ^at, Q, Ks^v are understood as the ones 
in the chiral limit. 

The Nb = 2 sector consist of terms with (local) two nucleon interactions, dibaryons and 
dibaryon- nucleon interactions. The terms with two nucleon interactions can be removed 
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by local field redefinitions [16, 17, 14] and will not be further considered. The terms with 
dibaryon fields and no nucleons in the rest frame of the dibaryons read 



i^D = ^0{p) + >C(9(i92) (2.5) 
where Co^p) is the 0{p) Lagrangian^ , 

'Dl[-id, + 5'^)lD^ +Dl(^-ido + S^^^^D, + zc,,(^DlTr{u^^^ (2.6) 



where Dg — D^Ta and 5^., i — s^v are the dibaryon residual masses, which must be much 
smaller than A^, otherwise the dibaryon should have been integrated out as the remaining 
resonances have. The negative signs of the time derivatives are chosen this way in order to 
eventually reproduce the signs of the effective range parameters. As discussed in [8], they 
do not imply any violation of unitarity. 

The covariant derivative for the scalar (isovector) dibaryon field is defined as doDg = 
doDs + l[[u,dou\,Ds]. Co{p2^ is the 0{p^) Lagrangian, 

Co(p2) ^siTr[Ds{uM^u + u^Mu^)DI] + S2Tr[Dl{uM^u + u^Mu^)Ds] + 

+ viDl . D^Tr[u^Mu^ + uM^u] + • • • ^^'^"^ 

5^, z = 1, 2, and Vi^ are low energy constants (LEG). We have only displayed here the terms 
which will eventually contribute in our calculations. The complete list of operators is given 
in the Appendix B. 

The dibaryon- nucleon interactions will also be needed at NLO, 

^DN — ^D^N + ^DTV + • • • (2-8) 

At LO they read 



4jv =^{N^aWN*)D,^+^{N^a'r^r^N)Dl+ 
V 2 v 2 



(2.9) 



— 1/2 

Ai ' ^ i — s^v^ and at NLO 



+ ^(7Vt^2^^2^2^*^ . ^ ^{N^r\''aD''N) ■ 31 (2.10) 

V 2 V 2 

+ ^{DiN^r''a'a^DjN*)Di + ^{DiN^T''a'^a'DjN)Di^ 

V 2 v 2 



o — ao 



Figure 1: The dibaryon propagator gets an important contribution from resuming the bubble self-energy 
diagrams 

Again, we have only displayed here the terms which will eventually contribute in our calcu- 
lations. The complete list of operators is given in the Appendix B. 

As discussed in [8], the dibaryon field propagator gets an important contribution to the 
self-energy due to the interaction with the nucleons (Fig.l), which is always parametrically 
larger than the energy E. As a consequence the LO expression for the dibaryon field propa- 
gator becomes (in dimensional regularization (DR) and minimal subtraction (MS) scheme), 

J-s,v, (2.11) 



{p — ^/ttinE) rather than the tree level expression i/{—E + 5'^. — iif). The size of the 
residual mass can be extracted computing the LO amplitude using the propagator (2.11) 
and matching the result to the effective range expansion (ERE), 

2 

ai are the scattering lengths. Therefore for p ^ the full propagator can be expanded. 
The first term of this expansion will be the LO propagator (Fig. 2), 

^ i = s,v, (2.13) 



AjrriNP 

the second term will be an effective vertex taking into account the effects due to i(— £^ + 5^.). 
Higher order terms in this expansion will be equivalent to multiple insertions of this vertex. 

Furthermore for p ^ 5^. the LO Lagrangian becomes both scale and SU{4:) (spin- 
fiavor Wigner symmetric) invariant, if the interactions with pions are neglected [18]. Indeed, 
concerning 5'C/(4), the single nucleon sector is obviously invariant. Moreover, since all terms 
in (2.6) become subleading, one can redefine the dibaryon fields in such a way that all 
couplings in (2.9) are equal. In that case the dibaryon-nucleon interactions become SU{4:) 
invariant if the two dibaryon fields are chosen to form a 6* representation of SU{4:). Scale 
invariance also holds because the dibaryon fields only appear in (2.9) and their scaling 
transformations can be chosen such that those terms are invariant. 



^The last term starts contributing to the nucleon- nucleon amplitudes at N^LO and, hence, was not 
displayed in [8] 
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Figure 2: The LO dibaryon propagator for p Srm is the first term in the expansion of the fuh dibaryon 
propagator around {—E + 5^.) = 0. The second term is an effective vertex. 

Moreover equation (2.11) implies that the dibaryon field should not be integrated out 
unless p <C instead of <C S^^. as the tree level expression suggests. If 5^^. <C m^^ it 
should also be kept as an explicit degree of freedom in the ^NNEFT, like in Refs. [19, 20, 21]. 
Nevertheless, if the dibaryon fields are integrated out, one can still organize the calculation 
in terms of nucleon fields by taking into acount suitable correlated enhancements in the local 
four nucleon interactions [13]. This is due to the fact that the path integral over dibaryon 
fields is Gaussian and can be carried out exactly. 

Except for the above mentioned contributions to the self-energy of the dibaryon fields, 
which become LO, the calculation can be organized perturbatively in powers of 1/A^. Hence 
one expects that any UV divergence arising in higher order calculations will be absorbed in 
a low energy constant of a higher dimensional operator built out of nucleon, dibaryon and 
pion fields (note that the linear divergence in the self-energy of the dibaryon fields due to 
the diagram in Fig. lb can be absorbed in S^.). 

We shall restrict ourselves in the following to energies E < m^/A^ <C ttIt^^ which implies 
nucleon three momenta ^ 171^^. We shall follow the strategy of [22], which was inspired in 
the formalism of [23], and shall build a lower energy EFT, pNNEFT, with no explicit pion 
fields: the effects due to the pions will be encoded in the potentials (and redefinitions of the 



3. The potential nucleon- nucleon effective theory 
with dibaryon fields 

For energies E ^ m^/A^ <C m^^ the pion fields can be integrated out. This integration 
produces nucleon- nucleon potentials and redefinitions of low energy constants. Since 5^. ^ 
m^/A^ the dibaryon fields must be kept as explicit degrees of freedom in pNNEFT. 
The Lagrangian in the A^^^ = 1 sector reads 



In the Nb = 2 sector further two nucleon interactions (potentials) are induced. They 



LECs). 




(3.1) 
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read, 

Cnn =\ j dhN^a"r''N{x,)Va0,,,ixi - X2)N^ a^r'^ N{x2) (3.2) 

a:^ = = r = xi — X2 and x = {xi + X2)/2 where Va/s-paixi — 2:2) is a generic potential 
(a, /3, p, cr = 0,1,2,3; t° = (7° = 1), which may be calculated in an expansion in 1/A^ (in 
fact, beyond one loop it becomes an expansion in [24]) 
The terms with dibaryon fields and no nucleons read, 

J^'d = Ka {-^do + 5^}jD1 + Dl {-ido + 5^^)d, (3.3) 

5^., i — s^v are the (redefined) dibaryon residual masses. Note that Srrijsf in (3.1) can be 
reshuffled into 5^. by local field redefinitions. Note also that because of 5^. <C m^^ the quark 
mass dependence of 5 mi is a leading order effect. 

The dibaryon- nucleon interactions remain the same as in (2.8), except for the values of 
the Ai which get modified. 

A ^ A ^ (3.4) 

+ ^^(TVtT^aa^iV*) • + ^(TV^rV^a^TV) . Dl 
V 2 V 2 

+ ^(NW^^a^d'^N*) • + ^{N^T'^a^ad'^N) • Dl (3.5) 

V 2 V 2 

+ ^{diN^r^a'a^djN*)Di + ^{diN^T^a^a'djN)Dl^ 
v 2 V 2 

The calculations in pNNEFT can be organized in ratios E/p and p/A^ (recall m^j- ^ p). 

The UV divergences arising at higher orders will be absorbed by local terms build out of 

nucleon and dibaryon fields. 



4. Matching NNEFT to pNNEFT 

In the Nb = 1 sector one loop pion contributions produce energy independent terms which 
are 0(m^/A^) [25] and hence relevant for the N^LO calculation, which together with the 
contribution 0{m^/A^) from terms proportional to the quark masses make up the nucleon 
residual mass Sttin in (3.1). 

In the Nb = 2 sector the dibaryon residual masses also get contributions 0{m^/A^) 
from terms proportional to the quark masses in (2.7), and 0(m^/A^) ones from higher loop 
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(a) (b) 

Fi gure 3: Order 0(m^^^/A^^^) contributions to the dibaryon residual mass with one radiation pion. (a) 
These three diagrams sum up zero by Wigner symmetry, (b) Wigner symmetry is violated by (2.6), this is 
by insertions of i{—E + Smi)- Naively we would expect these diagrams to be of higher order, 0(m^^^/A^^^), 
but the energy term is enhanced by the radiation pion up to 0(m^/A^). Hence the cross in this diagram 
stands only for an insertion of the energy. 




Figure 4: Order 0(m^/A^) contributions to the dibaryon residual Figure 5: Order 0(m^/ A^) con- 
mass with one radiation pion and one potential pion. Only diagrams tribution to the dibaryon residual 
with the potential pion inside the radiation pion loop contribute [10]. mass involving the Csv vertex 



diagrams involving radiation pions, like the ones in fig. 4, fig. 3b and fig. 5, 
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c 



s 



2/27 \ 4n J ' \2f^J87rAl '''Kf^JSnA, 



2PJ V 47r y ■ \2PJMAI '"^K^J Stt A, 

(4.1) 

One may have also expected a contribution 0(m^^^/A^^^) from diagrams in fig. 3a, but 
they add up to zero. This is not accidental, but a non-trivial consequence of Wigner sym- 
metry [18]. 

The dibaryon-nucleon vertices may in principle get 0(m^/A^) from a pion loop, but they 
turn out to vanish, except for those which reduce to iterations of the OPE potential which 
will already be included in the calculations in pNNEFT and must not be considered in the 
matching. Note that for this to be so the matching calculation must be done according to 
the prescriptions of ref. [26], which we briefiy recall in the Appendix A. This prescription 
gives results which differ from the on-shell prescription of ref. [10] and are usually simpler. 
Agreement is eventually recovered at the level of physical amplitudes, in which a number of 
cancellations occur for the on-shell prescription. There is, however, a two loop contribution 
of this order involving radiation pions from the diagram in fig. 6. 
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Figure 6: Matching of the effective vertex of the pNNEFT theory with the NNEFT vertex diagram. 




Figure 7: LO diagram 

Finally, in the two nucleon interactions (3.2), the one pion exchange is the only rele- 
vant contribution at this order, which produces the well known one pion exchange (OPE) 
potential, 

- = -|| / ^^M^<5-e--(---) (4.3) 
2/2 J (27r)3 q^ + ml 

for a, P, p,a = 1, 2, 3 and zero otherwise. 



5. Calculation in pNNEFT 

When p ^ mj^^ we have already integrated out all higher energy and momentum scales in 
pNNEFT, and hence we already have the optimal EFT to carry out calculations. Moreover, 
for this momentum both the time derivative and the residual mass in the dibaryon Lagrangian 
are small and can be treated as 0{m^/A^) perturbations. 

Let us then focus on the calculation of nucleon- nucleon amplitudes up to N^LO. At LO 
we get from fig. 7 the following S wave scale covariant Wigner symmetric amplitudes 

4:7T 

^io = ^ j = s,v (5.1) 

ttinP 



At NLO we get from the diagrams in fig. 8a, 

"8/2 V 47r / ' m2 



^NLO = -^[ H^ + Z::^l(^Lo) i = s,v (5.2) 
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V 




(a) 



(b) 



Figure 8: (a) NLO diagrams with one potential pion exchange, (b) NLO diagram with one i{—E + ^^.) 
insertion 




Figure 9: N^LO diagrams with two potential pion exchange 



and from the diagram in fig. 8b 



J: 



NLO 



4AJ 



i^lo) i = s,v 



At N^LO we obtain the following contributions. Prom the diagrams in fig. 9, 



(5.3) 



mt 



2 V Vm, 



2ip 



+ 



llm(L.J "^- + '^M ) + llnfl + arctanf^) 
4 V \-m^ + 2ipJ J 2 V mt/ \m^J 



(5.4) 




Figure 10: Diagrams with one vertex from (2.10). 
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(5.5) 

The sum of the diagrams in fig. 10 turns out to be zero for the ^Sq and channels (they 
only contribute to the ^^i-^Di mixing, see below). This can be understood as follows: these 
diagrams involve corrections to the nucleon-dibaryon vertices of order (^)^. We can redefine 
the dibaryon fields in order to remove these corrections from nucleon-dibaryon vertices, as a 
consequence these corrections would appear in the Nb = 2 sector (3.1), however, since this 
sector is subleading, the new operators induced by the field redefinition in this sector are of 
higher order. 

From the diagrams in fig. 11 




(5.6) 

Finally there are two contributions coming from relativistic corrections. The first one 
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Figure 12: The cross in one of tlie nucleon propagators stands for tiie use of tiie relativistic correction. In 
order to compute this diagram the relativistic correction has been expanded up to first order 

comes from using 2^"^ s^"^^^) instead ofi/{p^ — ^^+ie) as the nucleon propagator. 

We obtain the contribution in Fig. 12, 

Another contribution arise when using the first relativistic correction to the dispersion 
relation of the nucleons, 

P' = 7^- (5-8) 
^ 2mN 8m%' ^ ' 

in the bubble self-energy diagram of Fig.l. This results in the following contribution to the 
amplitude, 

However a new normalization of the amplitude that takes into account the new dispersion 
relation has to be considered, 



3 

this new normahzation induces a new contribution to the S matrix, i ^^^^ ^-^^, that exactly 
cancels the contributions coming from the aforementioned relativistic contributions to the 
amplitude. 

The ^^i-^Di mixing amplitude has no contribution at LO. At NLO the first two diagrams 
of Fig.Sa are the only contribution. 
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At N^LO diagrams of Fig. 9 with one (or two) potential pion exchange in the nucleon 
external legs give the following contribution, 



(5.12) 



The X, y, Z functions where defined in [11] and we write them here for completeness, 



X{a) 
y{cl) 



3 ?)% ia . / 1 3 



ln(l - 2ia) , 



2 

~5 ^ 10^2 



3 / 3 5 



2a 

8a^ ' 4a^ 5 



tan-^(a)- (A + _l_)tan-i(2a) 



(15-4a2)^ 2^ (3 + 16a2 + 16a4)^ 
+- — T^TT-^ — - ln(l + q; ) - ^ T:r-^ ^Im 



+i 



80tt6 

3 J_ _ « _ (3 + lOa^) 
8tt3 + 2a~ 2 16tt5 



8a^ 4q;3 
20:^ + 



(5.13) 
(5.14) 



' I l + 4a2 ) + Lt2{-2a - za) 



ln(l + 4q;2) + — — ln\l + Aa^) 



Z{a) 



128^7 
9 15 



7 9i 21 3i 3ia 29a'^ /3a^ 
~40 ^ T6tt3 + 40tt2 ^ iOtt ~ ~ ^ ^OO" ^ ~ T6a4 ~ So^ 
3 (16«^ - 50tt3 - 4z«2 - 15a + 15i) , , 

^(-9i + 27«-2W + 78c»-16«') ^^^^ _ ^^^^ 



In 2 



(5.15) 



(9 + 48^2 ^ 48^4^ 



- In' 1 - 2ia) + 2Lz2 -1 + 2ia) + Li^ , , ^. + — 
2 V — 1 + ZiaJ 4 



The derivative vertex of (2.10) proportional to B'^ contributes to the mixing amplitude 
through the first diagram of Fig. 10, 



TDI 

AmixJI j./Or?^ AV.^ 



(5.16) 



The last contribution to the ^Si-^Di mixing amplitude comes from the first diagram of 
Figlla, 



A mix, III ./7^( 9a\ 



g\ \mNm^ /-E + S^^ 



VP 47r 



4A2 



3/m^\2 /3/m^\3 ( 2p 



<-)■ 



+ i 



3m„ 



+ 



P 



+ 



4p 2m^ V 4p 16 V 



(5.17) 
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The ^Di amplitude starts at NLO with the contribution coming from the one pion exchange 
diagram, 




(5.18) 



At N^LO there are two contributions from Fig.9 from the two diagrams in which all 
external nucleon legs have a potential pion exchange. The corresponding amplitudes are, 



'Dij ( 9a Y { mNm^ \^ ( 
■N^LO [2f2) [ 4^ 




6. The pionless nucleon- nucleon effective field theory 

2 

For P ^ ^ the calculation must be organized in a different way. This is very much facilitated 
if we integrate out nucleon three momenta of the order of first, which leads to the so called 
pionless nucleon-nucleon EFT (/NNEFT) [12, 13, 14]. This EFT has been successfully used 
in numerous processes at very low energy (see [7] for a recent review). The Lagrangian of the 
Nb = 1 sector of this theory remains the same as in pNNEFT (3.1) (the relativistic correction 
becomes negligible) . For the A^^ = 2 sector the only formal difference from pNNEFT is that 
the non-local potentials (3.2) become local and can be organized in powers of /m^. The 
OPE potential in (4.3) becomes 0{p^ /rin?^K^ and hence beyond N^LO in this region. The 
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derivative dibaryon-nucleon vertices in (3.5) also become beyond this order. The remaining 
terms in the Lagrangian are the same as those in pNNEFT, namely (3.3) and (3.4), with 
the parameters redefined as follows. Diagrams in Fig. 8a and Fig. 9 containing one (or two) 
potential pion inside a nucleon bubble will contribute to the dibaryon time derivative term 
as well as the dibaryon residual mass. Contributions to the dibaryon time derivative can 
be reabsorbed by field redefinitions of dibaryon fields, while contributions to the residual 
mass simply redefine it. The dibaryon-nucleon vertex (3.4) gets contributions from diagrams 
containing one (or two) potential pion in the dibaryon-nucleon vertex, redefining the Ai. 
There are also higher order self-energy diagrams for the dibaryon fields not shown in paper 
which contribute to the redefinitions of the residual dibaryon masses at the order we are 
interested in, like the ones involving three OPE or the two pion exchange potential in a 
nucleon bubble. We will chose to reshuffle all matching contributions to the dibaryon-nucleon 
vertices to the residual masses through field redefinitions of the dibaryon fields. This way the 
coupling constants will remain the same as in pNNEFT while all the new dependences 
are carried by the residual masses. 

Since the dibaryon residual masses are no longer small, but of the same order, when 
compared to residual masses have to be kept in the dibaryon propagators. Hence we will 
use (2.11) rather than (2.13) as the dibaryon propagator. 

The LO amplitude for the pionless EFT is obtained from the diagram in Fig. 7 using the 
new dibaryon propagator, 

A';* = -I j = 5, V. (6.1) 

dm + 

Note that the LO amplitude is of order 0(1/ ml.) instead of 0{l/m^A^) as in pNNEFT, 
however since contributions to the S matrix are proportional to the momentum the final size 
of the LO contributions to the observables (as well as the NLO and N^LO ones) remains 
the same as in the high energy region. Note also that both scale invariance and Wigner 
symmetry are lost in the low energy region. 

The form of the amplitude remains the same at NLO (only Ai and get redefined). At 
N^LO (i.e. 0(1/A^)) a contribution corresponding to Fig. 8b arises. 



K*^LO--(i^]i^for ^ = s,v. (6.2) 



E 
Ml 

The form of the N^LO expression turns out to be valid also up to N^LO (i.e. 0{m^/A^)^ 
again only Ai and get redefined) The sum of diagrams in Fig. 10. is no longer zero but 
the momentum dependence of the vertex involved makes them beyond N^LO. 
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Af^^(MeV-V2) 


^^(MeV-V2) 








^(MeV) 


NLO 


0.0291 




-1.40 




-3.90 




N^LO 


0.0361 


0.0277 


-13.7 


-17.4 


2.10 


-1.89 



Table 1: Fit values of the parameters for the ^S'o channel. 



No contributions to or to jt^^'^ appear up to N^LO (the first diagram of Fig. 11. a 

contributes to at N^LO; this amplitude matches a straightforward expansion for 

p <C ttIt^ of the pNNEFT mixing amplitude) . 

7. The ^S{) channel 

In this section we compare our results for the ^S'o channel with its corresponding phase shift 
data. In order to compute the phase shift the amplitude has been introduced in exp(2i5) = 
1 + ipm]s[A/27r. After expanding both sides in powers of (rrij^/A^)'^ the expressions for 5^^, 
5^^^ and, 5^'^^ are obtained. 

We will not display the results for 6^^. At this order our approach does not uniquely 
determine the phase shift in the high energy region. This can be easily seen if the expression 
for the phase shift is expressed in terms of the real and imaginary parts of the amplitude, 

Since our LO amplitude (5.1) has no real part, then 6 — ±7r/2. Continuity with the low 
energy expression selects the plus sign. 

As and 6ms receive corrections in the matching from NNEFT to pNNEFT, both at NLO 
and N^LO. If the whole expressions for Ag and were to be used in the N^LO amplitude, 
higher order terms would be introduced. Therefore we will differentiate between A^^^ and 
Af^^ as well as between 6^^^ and 5^'^^, which we will consider as independent parameters. 
Recall that the expression for the phase shift in the low energy region shares the same Ag 
as in the high energy one, but has an independent 6ms •> which we will label 6ms • Because of 
this shared parameter [Ag) we have made a common fit of the low and high energy region 
phase shift at each order. The low energy region phase shift (calculated in ^fNNEFT) has 
been fitted to data in the 0-3MeV range and the high energy region phase shift (calculated in 
pNNEFT) to data in the 3-50MeV. Results for the channel parameters are summarized 
in Table 1. An alternative fitting procedure was presented in [27]. 

The phase shifts are plotted in fig.l3(NL0) and fig.l4(N^L0) versus center of mass (CM) 
energy. The low energy region and the high energy region phase shifts have been plotted in 
the 0-4MeV and l-50MeV range respectively. 
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Figure 13: '(Color online) 'Plot of the NLO expression for the ^S'o phase shift versus CM energy. The 
blue Hue shows the Nijmegen data for the ^5*0 phase shift, while the red and green line correspond to the 
high energy and low energy expressions respectively. The fitting procedure is explained in the text. Error 
bands correspond to for p < and to for p > 
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Figure 14: '(Color online) 'Plot of the N^LO expression for the ^5'o phase shift versus CM energy. As in 
the previous figure the blue line shows the Nijmegen data for the ^5*0 phase shift. The red line corresponds 
to the high energy expression and the green one to the low energy one (which totally overlaps the data). 
Error bands correspond to for p < and to for p > 
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C^'^(MeV) 


Sm^ '^(MeV) 


BUA, (MeV-2) 


NLO 


0.0305 




12.14 




8.30 






N^LO 


0.0431 


0.0429 


-9.29 


-13.3 


-23.0 


19.9 


-1.78 • 10-^ 



Table 2: Fit values of the parameters for the — ^ Di channel, excluding the N^LO ^S'l phase shift in 
the high energy region. 



8. The ^Si-^Di channel 

In this section we analyze ^Si-^Di channel. We compare the ^Si and ^Di phase shifts to 
data as well as the mixing angle. The usual expression for the S-matrix in this channel, 



.prriN ( ^'"^^ \ _ / e^^^*'^^' cos(2e) ze'^*'^^'+^^*'°^) sin(2e) 

^ 1^ A'^^ )~[ ie^5<^^^'+^5<^^^' sin(2e) e^^^*'"^' cos(26) 

(8.1) 

To obtain the phase shift expression at each order we expand both sides in powers 
(m^/A^)^^ as we did in the previous section, and solve the resulting system to obtain S^'^^ ^ 
5-'™ and, 5-'^'^^; and, ^^^i^^^^^; e^^^ ^n-lo^ ^here is no e^^ or 5'^^'^^ 

due to the fact that ^^'^ and J^^^ start at NLO. 



The fitting procedure for the NLO result is analogous to the one used for the ^S'o channel. 
For the N^LO one, several changes had to be introduced. A common fit to the low energy 
phase shift and to the mixing angle have been made, whereas the high energy phase shift has 
been left out and fitted independently. This is because all attempts to fit the high energy 
phase shift together with the other two expressions failed. The N^LO pNNEFT phase shift 
fit delivers a value for A^^^^ (Table 3) which is far away from the expected natural size. We 
think this is the reason why we were unable to perform a successful common fit: whereas the 
mixing angle and the low energy phase shift favor natural size parameters, the high energy 
phase shift does not. This is a clear sign that our approach fails to converge in the ^Si-^Di 
channel, we will comment on it further in the next section. Note that the parameter B'^/A^ 
only appears in the N^LO the mixing angle, e^^^, S^^^'^^^ and, S^^i^^'^^^ do not contain 
free parameters. An alternative fitting procedure was presented in [27]. 

The ^S'l phase shifts are plotted in fig.l5(NL0) and fig.l6(N^L0), the mixing angle in 
fig. 17 and the ^Di phase shift in fig. 18. All ^Si-^Di channel plots are versus CM energy. The 
low energy region and high energy region phase shifts have been plotted in the 0-4MeV and 
l-50MeV range respectively, the mixing angle and the ^Di phase shift have been plotted in 
the 0-50MeV range. Results for the ^Si-^Di channel parameters are summarized in Table 2 
and Table 3. 
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Figure 15: '(Color online) 'Plot of the NLO expression for the ^Si phase shift versus CM energy. The 
blue line shows the Nijmegen data for the phase shift, the red line corresponds to the high energy region 
expression and the green to the low energy region one. The fitting procedure is explained in the text. Error 
bands correspond to for p < and to for p > 
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Figure 16: '(Color online) 'Plot of the N^LO expression for the phase shift versus CM energy. The 
blue curve is the Nijmegen data for the ^Si phase shift, while red line corresponds to the high energy region 
expression. The curve for the low energy expression totally overlaps with data. The fitting procedure is 
explained in the text. Error bands correspond to for p < m^^ and to for p > m^^ 







,5^fO(MeV) 




0.0206 


0.00996 


35.3 


3.04 



Table 3: Fit values of the parameters delivered by the N^LO phase shift in the high energy region 



18 



e 

^(deg) 
10 

8 

6 

4 

2 



-2 

Figure 17: '(Color online) 'Plot of the mixing angle versus CM energy. The blue line shows the Nijmegen 
data, the green and red lines the NLO and N^LO expression respectively. The NLO expression has no free 
parameters. The free parameters of the N^LO expression have been fitted as explained in the text. The light 
green (light red) error bands correspond to ±(^^)^^^^ for p < and to ±(^^)^(^^ for p > m^. 
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Figure 18: '(Color online) 'Plot of the ^Di phase shift versus CM energy. The blue line shows the 
Nijmegen data, the green and red lines the NLO and N^LO expression respectively. Neither the NLO or the 
N^LO expression have free parameters. The light green (light red) error bands correspond to ±(^^)^*^^'' for 
p < and to ±(^)^(^) for p > m^. 
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9. Conclusions 



We have calculated the nucleon-nucleon scattering amplitudes for energies smaller than 
the pion mass in the ^Sq and the ^Si —^Di channels at N^LO in a chiral effective field theory 
which contains dibaryon fields as fundamental degrees of freedom, the NNEFT. The large 
scattering lengths in the ^S'o and the ^S'l channels force the dibaryon residual masses to be 
much smaller than the pion mass. We organize the calculation in a series of effective theories, 
which are obtained by sequentially integrating out higher energy and momentum scales. We 
first integrate out energy scales of the order of the pion mass. This leads to an effective 
theory with dibaryon and nucleon fields, pNNEFT. The latter interact through potentials. 
For three momenta of the order of the pion mass, the scattering amplitudes are calculated in 
this effective theory. For three momenta much smaller than the pion mass, it is convenient to 
further integrate out three momenta of the order of pion mass, which leads to the ^-NNEFT, 
and carry out the calculations in the latter. By splitting the calculation in this way we 
can take advantage of the modern techniques of the threshold expansions and dimensional 
regularization so that all integrals only depend on a single scale [28, 29, 26]. There is no 
need to introduce a PDS scheme [9]. The technical complexity of the N^LO calculation is 
similar to the one in the KSW scheme [10], but our final expressions are simpler. 

The numerical results for the phase shifts and mixing angle are also similar to the ones 
obtained in the KSW approach. Hence a good description of the ^Sq channel is obtained 
, but for the ^Si Di channel our results also fail to describe data. The ^Si phase shift 
shows a good agreement with data up to center of mass energies of bOMeV at NLO, but at 
N^LO the range of the agreement is reduced, up to 20MeV. only, even when the high energy 
region of this channel is fitted independently, as in the plot of fig. 16. The mixing angle 
poorly agrees with data, but shows a marginal improvement from NLO to N^LO. The N^LO 
mixing angle plot is significantly different from the one shown in [10], this is a consequence 
of making a common fit of the ^S'l phase shift in the low energy region and the mixing angle. 
For a different fitting approach with plot closer to [10] see [27]. Comparison with data for 
the ^Di phase shift it is never good. Particularly worrying is the fact that for the ^Si and 
the ^Di phase shift the N^LO calculation compares worse to data than the NLO one. The 
reasons of this failure can be traced back to the iteration of the OPE potential , the first 
diagram in Fig.9, which gives a very large contribution [10]. This may be interpreted as 
an indication that pion exchanges must be iterated at all orders, as originally proposed by 
Weinberg [1]. However, the removal of the cut-off in this approach appears to require an 
infinite number of counterterms, one for each partial wave [30, 31, 32] (see also [33]). A 
very recent proposal, which keeps the essentials of KSW counting, consist in introducing a 
Pauli-Villars regularization for the pion exchanges and staying at the regularized level [34]. 
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This seems to produce slightly better results, but it is unclear at the moment that, staying at 
regularized level, this approach is superior to Weinberg's one [7] (see [35, 36] for very recent 
efforts on the renormalization of Weinberg's approach). 

Before closing we would like to add a few remarks to the current situation. Starting 
at N^LO the expansion parameter in our approach is ^/pJX^ rather than p/A^ (in fact 
it is an accident due to Wigner symmetry that up to N^LO the expansion parameter is 
the latter). The N^LO calculation, i.e. 0{vrij^ j would be relatively simple in our 
approach. Since fractional powers only arise from diagrams involving radiation pions, it 
would basically consist of taking into account E/ttIt^ corrections to the diagrams in Fig.3.a, 
and considering a further potential pion exchange or and (extra) internal energy insertion (a 
cross) in the diagrams in Fig.S.b and Fig.4. The outcome of the calculation, however, would 
amount to redefinitions of previously existing parameters, and hence it would not produce 
an improvement in the description of data. Qualitative changes are expected at N^LO. It is 
at this order, for instance, that the two pion exchange potential first enters the calculation. 
The N^LO calculation appears feasible in our formalism, but would require a major effort. 

Acknowledgments 

We thank Enrique Ruiz Arriola for explanations on refs. [30, 32]. We have benefited from 
the participation in the ECT* workshop "Bound States and Resonances in Effective Field 
Theories. We acknowledge financial support from the CIRIT grant 2005SGR00564 (Cat- 
alonia), MICINN grants FPA2007-60275/, FPA2007-66665-C02-01/ (Spain), the RTN Fla- 
vianet MRTN-CT-2006-035482 and the European Community-Research Infrastructure Inte- 
grating Activity "Study of Strongly Interacting Matter" (acronym HadronPhysics2, Grant 
Agreement n. 227431) (EU). JT has been supported by a FI grant from Departament 
d'Universitats, Recerca i Societat de la Informacio of the Generalitat de Catalunya and a 
MICINN FPU grant ref.AP2007-01002. 

A. Matching prescription 

In the matching calculations between NNEFT and pNNEFT there are regions in the integrals 
in which fc^ ~ k ^ ttIt^. In these regions the kinetic term of the nucleons is parametrically 
smaller than the energy and hence, following the ideas of the threshold expansions [28], it 
must be expanded. However, when one does so in the two nucleon sector one often encounters 
pinch singularities. It was argued in [26] that the pinch singularities can be rearranged in 
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such a form that they exactly match the pinch singularities of the effective theory when 
the kinetic term is also expanded in the latter. Consequently, once this rearrangement is 
done, pinch singularities can be safely ignored. As an example, let us consider the one loop 
contribution to the nucleon-dibaryon vertices produced by a pion exchange. External three 
momenta are of the order of m^^ and external energies much smaller than it. Once small 
scales are expanded we are faced with the following integral, 

r d^k k'y 1 1 _ _ r d^k k'y i i 

J (27r)^ A;2 - m2 + ie P + ie-k^ + ie~~ J (27r)^ + P + ie -P + ie ^ ' ^ 

The last expression matches exactly the contribution of the OPE potential to the vertex in 
pNNEFT, if the kinetic terms of the nucleons are correspondingly expanded. Hence, in this 
case there is no contribution to the matching. 

Note also that with this prescription the size of each diagram in the fc^ ^ k ^ 171^^ region 
is easily estimate since the integrals depend on a single scale. For instance, the first diagram 
in Fig. 3 (a) has also a contribution in this region but it can be easily seen to be higher 
order. This is not so if the on-shell prescription is used [10]. 



B. The complete NLO Lagrangian in the Nb = 2 sector 

We list here all the operators of the NLO Lagrangian in the A^^ = 2, many of which do 
not contribute to our calculations. We use for organization purposes the standard chiral 
counting, namely do di — 0{p) and the quark mass matrix ^A — 0{p^) 

2.1 The dibaryon Lagrangian 

The full list of operators in Co(p'^) of (2.7) follows 

Tr[Ds{uM^u + u^Mu^)Dl] , Tr[Dl{uM^u + u^Mu^)Ds] , 

Tr[DlDsUoUo] , Tr[DlDsUiUi] , Tr[DsDluiUi] , Tr[DluoDsUo] , Tr[DluiDsUi] 

31 . D^Triu^Mu"^ + uM^u] , 31 • 3^Tr[uoUo] , 31 • 3^Tr[uiUi] 

{Dl^Di + DlDi^)Tr[u,u,] , Tr[Dlu x u]3, + h.c. (B.l) 

d3lTr[uoDs] + h.c. , 3lTr[uodDs] + h.c. 

Tr[dDldD,] , (d3l) (d3,^ , 3ld'3, 

In Ref. [8] terms mixing the scalar and vector dibaryon as well as terms with space derivatives 
on the dibaryon field were not displayed^. 

^The two additional terms that were displayed in [8] turn out to be redundant. 
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2.2 The dibaryon-nucleon vertex 

The full list of operators in of (2.10) follows (hermitian conjugates are omited) 



(N^a^D.T^N*) {Tr{uoUo) , Tr{uiUi) , Tr{uM^u + Mu^)) (B.2) 
[uqDsUq , UiDgUi , DguM^u , DgU^Mu' , u'Mu^Ds , uM^uDg) r^a'^N* 

N^T^a'a^N* {S'^Tr{uoUo) , S'^Tr{ukUk) , S'^Tr{uM^u + u^Mu'') , Tr(MiMj)) Di 
N^a'{u'Ds , , e'^%UkDs , e'^'^UjDsUk , e'^'^DsUjUk) r^a^N* 
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